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Abstract 

Entanglement evolution of two independent Jaynes-Cummings atoms without rotating-wave ap- 
proximation (RWA) is studied by an numerically exact approach. The previous results in the RWA 
are essentially modified in the strong coupling regime (g > 0.1), which has been reached in the 
recent experiments on the flux qubit coupled to the LC resonator. For the initial Bell state with 
anti-correlated spins, the entanglement sudden death (ESD) is absent in the RWA, but does appear 
in the present numerical calculation without RWA. Aperiodic entanglement evolution in the strong 
coupling regime is observed. The strong atom-cavity coupling facilitates the ESD. The sign of 
detuning play a essential role in the entanglement evolution for strong coupling, which is irrelevant 
in the RWA. An analytical results based on an unitary transformation are also given, which could 
not modify the RWA picture essentially. It is suggested that the activation of the photons may be 
the origin of the ESD. The present theoretical results could be applied to artificial atoms realized 
in recent experiments. 

PACS numbers: 
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I. INTRODUCTION 



Quantum entanglement, one of the most striking consequences of nonlocal quantum corre- 
lation, is fundamental in quantum physics both for understanding the nonlocality of quantum 
mechanics^ and its role in quantum computations and communications^]. The entangle- 
ment would undergo decoherence due to the unavoidable interaction with the environment. 
As a result, an initially entangled two-qubit system becomes totally disentangled after evolv- 
ing for a finite time. This phenomena is called entanglement sudden death (ESD)[3| and 
has been recently demonstrated experimentally [4(. 

Many works are devoted to the ESD of the qubits coupled to an environment that results 
in irreversible loss, and the rotating-wave approximation (RWA) is made on the interaction 



of the qubits with the field 



0,1 



The effect of the counter-rotating terms on the ESD has 



been less studied, even for the Jaynes-Cummings (JC) model [6] where the qubit interacting 
only with the single-cavity mode. In fact, some investigations have also focused on the 
storing of entanglement in such a system [7j. So entanglement dynamics for two independent 
JC atoms is also of fundamental interest, and has been well studied only in the RWA^- 11]. 

On the other hand, the JC model is also closely related to condensed matter physics 
recently. It can be realized in some solid-state systems recently, such as one Josephson 



12| , the superconducting quantum 



3. 



and the most recently 



charge qubit coupling to an electromagnetic resonator 
interference device coupled with a nanomechanical resonator 
LC resonator magnetically coupled to a superconducting qubit |15[|. In conventional quantum 
optics, the coupling between the two-level "natural" atom and the single bosonic mode is 
quite weak, RWA has been usually employed. With the progress of the fabrication, the 



artificial atoms may interact very strongly with on-chip resonant circuits 12Hl5j , RWA can 
not describe well the strong coupling regime [if]]. Therefore, it is highly desirable to explore 
the entanglement dynamics for two independent JC atoms without RWA. 

However, due to the consideration of the counter-rotating terms, the photon number is 
not conserved, so the photonic Fock space has infinite dimensions, any solution without 



16 



21] has been 



RWA is highly nontrivial. In the recent years, several non-RWA approaches 
proposed in a few contexts. Especially, by using extended bosonic coherent states, three of 
the present authors and a collaborator have solved the Dicke model without RWA exactly 
in the numerical sense 17j . 
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In this paper, we employ a numerically exact technique to solve JC model without RWA 
by means of extended bosonic coherent states. The correlations among bosons are added 
step by step until further corrections will not change the results. All eigenfunctions and 
eigenvalues can be obtained exactly. Based on the exact solutions to the single JC model, 
we can study the entanglement evolution of two JC atoms easily. Analytical results without 
RWA based on an unitary transformation are also presented. The paper is organized as 
follows. In Sec. II, the numerically exact solution to the JC model is proposed in detail 
and analytical results in terms of unitary transformation are also provided. The numerical 
results and discussions are given in Sec. III. The brief summary is presented finally in the 
last section. 



II. MODEL HAMILTON! AN 

The JC model can be written in the basis ^ |f) , ||) j where the first index denoting 
spin-up and the second denoting spin-down 

Hjc = yO" 2 + wa + « + A(a + a + )a x (1) 

where a + and a are the bosonic annihilation and creation operators of the cavity, A and 
uj are the frequencies of the atom and cavity, g is the atom-cavity coupling constant, and 
<7k(k = x, y, z) is the Pauli matrix of the two-level atoms. Here we set h — 1. The detuning 
is defined as 5 = ui — A. 

In this paper, we describe two approaches to solve the single JC model without RWA. 
One is the numerically exact approach in terms of extended bosonic coherent states, the 
other is an analytical approach based on a unitary transformation. 

Numerically exact approach - In order to facilitate the calculation, we use a transformed 
Hamiltonian with a rotation around the y— axis by an angle 7r/4, so <j z — > cr x ,a x — > —cr z . 
This rotation can also be described formally by the following transformation, 

H' JC = VH JC V + = -^a x + ua + a + X(a + a + )a z , (2) 

where 
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By introducing the new operators 

A = a + g,B = a- g,g = (3) 

we can observe that the linear term for the bosonic operator is removed, and only the number 
operators A + A and B + B are left. Therefore the wavefunction can be expanded in terms of 
these new operators as 



For A operator, we have 



M'»=l" c " k H (4) 

En=0 °2n \n) 1 



A n , . (a + g/u) n ln , 
Vnl \/n\ 



The property of B operator is the same. 

The Hamiltonian (2) remains unchanged under the transformation level 1 — > level 2 and 
a + (or a) — > — a + (or —a). So we can set c 2 „ = ±(— l) n C\ n (or c\ n = ±(— l) n c 2n )• The final 
more concise wavefunction is supposed as 

v(*)> =1 . , ... , . I • ( 6 ) 



±En=o(-l) n c n |n) B 



The Schro dinger equation is 



A Ntr 



U(m - g 2 )c m ± — ^ D mnC n = E {±) c m) (7) 
n=0 

where 

flm „ = exp(-2^) g (-1) \ m _^_ m - (8) 

The /— th eigenfunction \(p'®) is then obtained, which can also be expressed in the original 
basis of Hamiltonian (1) as 

M"> = I 1! I = ^ + k' (,) > • (9) 



The wavefunction at anytime then reads 

M 



Ht)) = h(l) e M-iE t t) \<P {l) (t)) , M = 2(N tr + 1), (10) 
i=i 



where h"' is the coefficient to be determined. 

For the two independent JC atoms, the eigenfunction of the total system is given by 



where <fi is the eigenfunction the JC model i(= 1, 2) with 



(11) 



\<Pi) 



K2) 



(12) 



The eigenvalue of the total system is E — Ei © E 2 . The j — th wavefunction for the total 
system can then be explicitly expressed in the basis ( \\"\) , |f J,) , |||) , |^} 



|V> (j) ) 



mi 



,j = 1,2,4(^ + 1). 



(13) 



The wavefunction of the total system at anytime then reads 



2M 



m)) = E f U) eM-iEjt) |V 0) ) , M = 2(iV ir + 1), 
i=i 



(14) 



where f^' is the coefficient to be determined. 

The initial two Bell states with anti-correlated and correlated spins, which are denoted 
by the Bell state 1 and 2 respectively for convenience, in the original base are the following 
if we use the column matrix 

/ 0|00) \ / cos«|00) \ 



VBell 



cos a 1 00) 
sin a 1 00) 

Vo|oo) J 



^ (2) \ 

VBell J 



0|00) 
0|00) 
Vsina|00) / 



(15) 



where 1 00) denotes the photon number state in the two JC atoms, /™ is determined by 



using [-0(0) ) 



For two-qubit states, entanglement can be quantified by the concurrence 
calculated from the following reduced density matrix 

N 

P = Tr ph (m)) (m\) = E ; (fc) / ( °e-^-^n, 

k,l=l 



22|. It can be 



(16) 



where 4x4 matrix II is determined by r0^)- The 4 eigenvalues Aj of the matrix p in 
decreasing order give the entanglement 



= max[0, VAT- V^2 ~ Vh ~ VM- (17) 
Unitary transformation approach - In order to treat JC model without RWA analyti- 



cally, we perform a unitary transformation 18] on Hamiltonian ([I]) to eliminate the counter- 
rotating wave term 

H s = e s He' s , 
H = H + H h H = y<7 2 + ua + d, H x = g(a + + a)(cr + + a J), 

S = — (a + — a)(a + + a_), 
a; 

where £ is a parameter to be determined. The transformed Hamiltonian can be expanded 
in terms of g up to the g 2 term (higher terms are neglected), and we have 



H S = H + H? + H* + 0(g 3 ), 



fff = # a + [S, H ] = g(-£ - £ + l)(a + a_ + aa + ) + g(--£ - £ + l)(a + o+ + a<r_), (18) 

If we choose £ = u/(A + u), the counter- rotating wave term can be eliminated, then we 
have the following renormalized JC Hamiltonian in a RWA form 

H s _ _^£l (Jz J rUja + a j r g ef ^ a + (T _j ra(T+ ^ (20) 

2g 2 
(A + co) 

9eff = 9 fcrk) • (22) 

where a constant is removed. The effective detunning can be expressed in terms of original 
detunning 5 as 



For 5 = 0, we have 



5 eff = ^K ff = ^-S° eff , £„ = g. (24) 
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Note that even for zero detunning, we have finite effective detunning after the unitary 
transformation. What is more, the non-RWA results within the unitary transformation 
approach is expected to be very close to the RWA ones in the weak coupling regime, because 
the effective g is not changed, and only modification is the transition frequency which is 
reduced by a small amount proportional to g 2 . 

Follow the derivation for the RWA case in Ref. Q], if the initial state is the Bell state 
1, the concurrence for two identical JC atoms without RWA by the unitary transformation 
can be easily obtained 



C AB (t) = |sin2a| [l - AN 2 sin 2 (z/t/2)] 



where 



N 



1 + 



1 2 



(25) 



(26) 



2g 

, 1+1/(1-5/"") - 



V 




2g 2 (oo -5) 
(2u - 5) 2 



2g 



1 + 1/ (1 - S/u) 



For zero detunning, we have 



2a/1 



_9_ 

N ' 



(27) 



(28) 



This is to say, within unitary transformation approach, there is no ESD from the initial Bell 
state 1, similar to the RWA case. 

Similarly, we can get the entanglement evolution from the initial Bell state 2 

Cab (i) = [l - 4iV 2 sin 2 (z/t/2)] (|sin2a| - 8iV 2 sin 2 (z/i/2) cos 2 a) , (29) 

The ESD occurs for 

|tana| < 4iV 2 sin 2 (z/t/2). (30) 

Note form Eq. (128]) that the value of iVo is very close to the RWA one even for g — 1, so it 
is expected that the non-RWA results within unitary transformation approach is nearly the 
same as that in RWA in a very large coupling regime. Moreover, the unitary transformation 
approach could not provide essentially different results, because final renormalized RWA 
Hamiltonian is derived for all calculations. In our numerically exact studies, the RWA-type 
Hamiltonian is not necessary, so some new results may be obtained. 



FIG. 1: (Color online) The concurrence for atom-atom entanglement with the initial atomic state: 
(a) Bell state 1 for a = ir/4 and (b) Bell state 2 for a = it/ 12. G = 2g. The corresponding right 
figures are the bird's-eye view. 



III. RESULTS AND DISCUSSIONS 



We first study the evolution of the concurrence for atom-atom entanglement for two JC 
atoms without RWA when the states are initiated with the given Bell states (a is fixed) ( 
( !T5|) . Actually, the dynamics of entanglement in bipartite quantum systems is sensitive to 



initial conditions si, |23[ 

We consider the two same JC atoms for zero detuning. The concurrence for atom-atom 
entanglement as a function of time can be calculated with the use of Eqs. ( TT3|) and (|T4|) . 
As show in Fig. []Ja), for small G = 2g, the concurrence evolve as cos 2 (Gt/2), just follow 
that in RWA [8J for Bell state 1. For G above 0.5, entanglement can drop to zero, which can 
not happen in the two JC atoms in RWA for Bell state 1. The entanglement can be revival 
irregularly. JC model is just N = 1 Dicke model. We attribute this aperiodic entanglement 
evolution to the quantum chaos in finite Dicke model 24j. Without RWA, the system is 



integral, the emergency of the quantum chaos is impossible, so regular behavior is always 



observed js, 1 10 1 . 

It is interesting that the area of the ESD become larger with increasing G. At very large 
coupling constant, the revival of entanglement will not happen, in sharp contrast with those 
observed in the RWA. Without RWA the rule of the transfer of entanglement between the 



two-qubit subsystems derived in 



-11 



may not hold due to presence of the counter-rotating 
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FIG. 2: (Color online) The concurrence for atom-atom entanglement with the initial atomic Bell 
state 1 (upper panel) and the Bell state 2 (down panel). From left to right column, g = 0.05,0.1, 
and 0.3. The numerically exact results and those by the unitary transformation are denoted by 
the blue and red curves. 




FIG. 3: (Color online) The concurrence for atom-atom entanglement with the initial atomic state: 
(a) Bell state 1 for a = tt/4 and (b) Bell state 2 for a = ir/12. g\ = 2g2,G = g\ + g2- The 
corresponding right figures are the bird's-eye view. 

terms. Without RWA, we argue that there is no entanglement invariant because the photon 
number is not conserved. 

For Bell state 2, the ESD can happen even in the RWA^]. Without RWA, for small G, 
the entanglement dynamics show the same behavior. As shown in Fig. [U(b), the area of the 
ESD become larger with increasing G, and more wider than that in Bell state 1. 

As we know, for weak coupling between the atom and the single bosonic mode, the RWA 
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is a good approximation, the non-RWA treatment is not necessary. Our results for the en- 
tanglement dynamics also provide a evidence for this point. However, for strong coupling 
case, the entanglement dynamics demonstrate that the non-RWA should be considered es- 
sentially. As the new progress in the fabrication, some artificial atoms are just strongly 
coupled with the bosonic field. 

The non-RWA JC model can be also treated by the unitary transformation approach. 
Is this approach really good enough in the strong coupling regime? We now present the 
analytical results based on this approach for zero detunning according to Eqs. (1251) and (1291) . 
which are shown in Fig. |2J The corresponding numerically exact ones are also exhibited for 
comparison. It is obvious that the analytical results essentially deviate from the exact ones at 
g — 0.1. Recent experiments on the LC resonator coupled to a flux qubit demonstrated that 
the system operates in the ultra-strong coupling regime g — 0.1 [l^, which crosses the limit 
of validity for RWA in the JC model. Since the analysis based on a unitary transformation 
could not essentially change the RWA results, the present numerically exact results are more 
necessary. 

We next consider the two JC atoms with different atom-cavity coupling for zero detuning. 
Without loss of generality, we choose g\ = 2g 2 . The concurrence for atom- atom entanglement 
as a function of time can also be calculated with the use of Eqs. (fl3l) and (fl4|) . As shown 
in Fig. |3]^a), for small G, the concurrence evolve as cos 2 (Gt/2), following that in RWA js] 
for Bell state 1. As G increases, entanglement can drop to zero and is almost not recovered, 
quite different from that in the two same JC atoms. We argue that two JC atoms with 
different coupling strength suppress the atom- atom entanglement. This point is confirmed 
by the evolution of the initial Bell state 2 as indicated in Fig. [3(b). No matter whether 
RWA is taken into account or not, the area of the ESD in the evolution of the initial Bell 
state 2 is larger than that in initial Bell state 1. 

For more detail, we study the effect of the coupling strength on evolution of the con- 
currence for atom- atom entanglement. Without loss of generality, we only consider the two 
same JC atoms. As exhibited in Figs. H] and El for weak coupling g\ < 10~ 3 , no ESD is 
observed. For g\ > 10~ 3 , the ESD appear, and its area is enlarged with g. At the same 
time, the periodicity of the entanglement evolution is destroyed. 

One natural question is what the mechanism is for the ESD in two remote JC models. 
The photons should influence the atomic state considerably. We then calculate the average 
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FIG. 4: (Color online) Histogram of the concurrence for atom-atom entanglement with the initial 
Bell state 1. From left to right, g = 10" 4 , 10" 2 , 1CT 1 (top panel) and g = 0.25,0.5,1 (bottom 
panel). G = 2g. 
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FIG. 5: (Color online) Histogram of the concurrence for atom-atom entanglement with the initial 
Bell state 2. From left to right, g = 10" 4 , 10" 2 , 10 _1 (top panel) and g = 0.25,0.5,1 (bottom 
panel). G = 2g. 

photon number during the evolution. To show its role in the entanglement dynamics, we 
plot the average photon number together with entanglement for two limit coupling cases 
g = 0.001 and 1 in Fig. Owith initial Bell state 1 and Fig. [7] with initial Bell state 2. It is 
very interesting that the curves for the entanglement and the average photon number always 
show opposite behavior. It is highly suggested that the average photon number suppress 
the entanglement between two atoms considerably and sensitively for both weak and strong 
coupling. 
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FIG. 6: The concurrence for atom-atom entanglement and the averaged photonic numbers for two 
identical JC atoms with the initial Bell state 1 for (a) g = 0.01 and (b) g=l. G = 2g, a = tt/4. 
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FIG. 7: The concurrence for atom-atom entanglement and the averaged photonic numbers for two 
identical JC atoms with the initial Bell state 2 for (a) g = 0.01 and (b) g=l. G = 2g, a = n/12. 
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Finally, we turn to the effect of detuning on the entanglement evolution. It is found in 
Ref. jlO] that the atom-atom entanglement in RWA usually increases with ratio between 
detuning and coupling constant \8\ / g and is irrelevant with the sign of detuning. In other 
words, the entanglement decreases with g for given detuning. Without RWA, it is not that 
case. After the unitary transformation, we have the renormalized JC Hamiltonian Eq. ( 12 Op 
in a RWA type. So entanglement evolution depends on |5 e //|, which is related to both the 
magnitude and the sign of original detuning \8\ shown in Eq. (121]) . By Eqs. (122]) and (|23|) . 
the effective coupling constant can also be expressed as g e ff = 2g/ [1 + 1/ (1 — 8/lj)] if the 
counter rotating-wave terms are taken into account after the unitary transformation. It 
follows that g e ff decreases with \S\ for positive detuning and increases with \8\ for nega- 
tive detuning. So in the strong coupling regime where the counter rotating-wave terms is 
required, the effect of detuning is very sensitive to the sign of detuning. The above discus- 
sions should be also suited qualitatively to the numerical exact results. We calculate the 
entanglement evolution of the initial Bell states 1 and 2 as a function of detuning. The 
results for different coupling constants ranging from weak to strong coupling regime are list 
in Fig. [HJ The entanglement evolution is symmetric with detuning for the weak coupling 
and becomes asymmetric with the increase of the coupling constant. In the weak coupling 
regime, as shown in the left column, the entanglement increase with the value of detuning, 
and independent of the sign, similar to those in RWA. With the increase of the coupling, 
it is however observed that the entanglement decreases with the magnitude of the positive 
detuning, and increases with the magnitude of the minus detuning, in contrast to the case of 
RWA. In the strong coupling regime, the positive detuning stabilizes the entanglement and 
facilitates its revival, whereas the negative detuning suppresses the entanglement, facilitates 
its death, and reduces the period of entanglement revival. 



IV. CONCLUSIONS 



In this paper, based on the exact solution for the single JC atoms, we are able to calculate 
exactly the entanglement evolution of the two independent JC atoms without RWA. The 
results are essentially different from RWA cases in the strong coupling regime. The analytical 
results based on an unitary transformation are also given. It can modify the RWA results 
and could not provide an essentially different ones. Initiated from the Bell state 1, the RWA 
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FIG. 8: Effect of detuning on the entanglement evolution of the Bell initial state 1 with a = 
7r/4(top panel) and Bell initial state 2 with a = n/12 (bottom panel). From left to right column, 
g = 10" 4 , 0.02, 0.1. 5 = u - A, G = 2g. 



results show no ESD. The present numerically exact calculations for the non-RWA model 
show that the ESD could not be avoided, and the periodicity of entanglement evolution is 
destroyed by the presence of additional photons. We also suggest that the photons may 
suppress the entanglement and is just the origin of the ESD. The effect of the detunning on 
the entanglement evolution is also investigated. It is observed that the sign of detunning 
play a essential role in the strong coupling regime. The present theoretical prediction would 
be tested in a experimental study of ESD where the artificial atoms are made of circuit 
QED 12Ml5| if operating in the ultra-strong coupling regime. 
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